


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 


DSpace Repository 


Theses and Dissertations 


l. Thesis and Dissertation Collection, all items 


1965 


Optimal filtering of radar data. 


Simoes, Carlos P. 


Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/26670 


Downloaded from NPS Archive: Calhoun 


ША DUDLEY 


W II KNOX 


LIBRARY 
http://www.nps.edu/library 






Calhoun is the Naval Postgraduate School's public access digital repository for 

research materials and institutional publications created by the NPS community. 

Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholarty author. 


Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 


NPS ARCHIVE 


1965 
SIMOES, C. 





OPTIMAL: FILTERING OF RADAR DATA 


CARLOS P. SIMOES 


DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 








Е 

















OPTIMAL FILTERING 
OF 


RADAR DATA 


Хх * +} Хх Ж Ж 


Carlos P. Simoes 





OPTIMAL FILTERING 


OF RADAR DATA 


by 
Carlos P. Simdes 
7 


Lieutenant, Portuguese Navy 


Submitted in partial fulfillment of 
the requirements for the degree of 


MASTER OF SCIENCE 
IN 
ENGINEERING ELECTRONICS 


United States Naval Postgraduate School 
Monterey, California 


1965 





U. S. Naval Posteradunte School 


OPTIMAL FILTERING 
OF RADAR DATA 
by 


Carlos P. Simöes 


This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 
ENGINEERING ELECTRONICS 
from the 


United States Naval Postgraduate School 





Е 
——— — 





ABSTRACT 


An optimal method of filtering radar data is analyzed and experimented. 

First the theoretical bases of Kalman's approach to filtering 
problems is presented. Then the suitability of its conclusions to 
the filtering of radar data is studied. The form for the filter needed 
in such a case is derived and, finally, a simulation is made attempting 
to show the real usefulness of Kalman's filter. 

It is concluded that the method considered presents a very prom- 
ising set of features which may lead to advantageous practical use. 

The filtering process is based on the recurrent solution of matrix 
formulas. To accomplish a reasonably fast computation of the filtered 
values the solution was programmed in the FORTRAN language and run in 


a 1604 CDC computer. 
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|. Introduction 


We will be concerned in this work with the filtering of radar data. 

To our knowledge the solution of this specific problem, in a major- 
ity of cases at least, was never attempted with the help of a deep theo- 
retical foundation. So far the ''«,@ Tracker'' seems to have conveyed 
almost all the possibilities of practical processing of radar data. 

With the new approach to the general problem of filtering announced 
by Kalman [1] a new frame of practical processing stood up to the engineer. 

The "' ος, (2 Tracker'' has no inherent means or need of utilizing the 
known statistics of the targets to be tracked in any mathematical or even 
rational way. 

The Kalman approach starts out by quantifying the known parameters 
of the signal generating process. Broadly speaking, the class of prob- 
lems attacked by Kalman consists of filtering and predicting the output 
of Gauss-Markov random processes when its parameters and statistics are 
known. The method described by Kalman being particularly suited to 
machine computation was immediately explored in more than one field 
of possible applications. 

One branch of the work done with the Kalman filter occupied control 
engineers. Another one gave results in the tracking of artificial 
satellites. 

Rauch [2] and Magill [3] followed Kalman with the analysis of 
cases in which the process parameters are random variables. 

it will be shown that Magill's Technical Report is completely 
suitable to provide a theoretical basis upon which a radar data filter 


can be built with fairly good results. 








2. Statement of thc problem. 


The processing of data provided by the various types of radars 
can be readily implemented with the digital computer. 

The processing device must be able to allocate the echoes received 
to one of the various targets being tracked, to decide on the appearance 
of a new target, to give numerical significance to the distance, azimuth 
and other characteristics detected by the radar, to store sequentially 
the information being obtained and to provide all the data stored to 
any unit or entity allowed to obtain it. 

This data will be considered to be present and available for 
filtering purposes. 

The data being stored needs to be handled in such a way that the 
various tracks maintained are filtered as accurately as possible. It 
is also considered that a prediction will have to be made and trans- 
mitted to some other system. 

The characteristics of the data provided by the computer storage 
for an eventual solution of the filtering problem are, obviously, of 
great importance. 

Admitting that the allocation of echoes to targets was done 
previously with total confidence, we do not need to consider but the 
filtering of one target track. 

This way for each target we will have stored data with the 
following characteristics, 

a) Numerical values assigned to distance, azimuth and altitude. 

b) Times at which the position was detected. 

c) If relevant, data on the position, speed and course of the 


vehicle where the radar is mounted. 








Adding to this data collected in the machine storage all the char- 
acteristics of the target defining its dynamic model are also necessary 
to improve the tracking efficiency. We may specify at this point that 
it is absolutely necessary to have some ''a priori' degree of knowledge 
about the dynamics of possible targets when the first point of the track 
is treated in the filtering process. As will be seen, the filter will 
not start its recurrent computation unless previous information on the 
characteristics of the possible targets is available. 

We need then: 

d) A dynamic model capable of describing the motion of the tracked 
object. 

e) Some statistics of the inputs to that model. 

Finally, the source of the noise to be filtered out has to be 
considered. Noise will be introduced by the radar itself, by the 
processing prior to the filter action, and by any anomalies in the 
propagation of the electromagnetic waves. 

Let us consider then: 

f) Statistics of the error introduced by the equipment. 

With this data available the filtering process will aim at 
reproducing the real track of the target, eliminating errors added 
during the acquisition of data. 

We will see that this objective is not completely achieved but 
that the mean square error of the positions will be reduced to a 
minimum, with the help of the Kalman filter. 

In d) we referred to the need of a dynamic model capable of describ- 
ing the motion of the object being tracked. More precisely we can say 
that this model is going to include or use, in one way or another, al! 


the data gathered by us or available in storage. We may even say that 
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this dynamic model is central to the solution we are seeking. With 
this in view let us proceed to formalize the model and obtain a more 
rigorous concept of the problem presented to us. 
The signal or message and the noise precesses are given by, 
5064) = (VS + 2 (ём) ц (е) 
УС) = Т) 5(2) 
W(tw)- Pult) WÉ) + Du ry) Uo (tri) 
n(t) «^ (0 w(t), 


Where the various symbols have the following meanings: 


and, 


SCH is the present state vector of the message process, i.e., 
a column matrix formed by the target components of distance, velocity, 
acceleration, etc., on some coordinate axis (which we assume to be the 
X-axis from now on). In chapter three we will indicate how the other 
components can be treated by the same filter. 

2 (2) is the transition matrix of the message process, interpreted 
as the matrix by which we must multipy the present state to obtain the 
next, as will clearly be seen in the end of this chapter. This includes 
the knowledge of d). 

U: (t) is a vector that includes all the possible inputs to the 
message process, as wind gusts (which we will consider as being able 
to change instantly the position of the target), air speed and varia- 
tions of acceleration as caused by the pilot's action. As later will 
be seen, these inputs have to be white gaussian random variables in 
order that the theory can be applied. 

Ds (€) is the distribution matrix of the message process and 
determines in which way the inputs affect the various state vector 


components. 





Y (E) is the message output which may consist of one or more of 
the state vector components. 
/(έ»/) is the output vector of the message process which determines 
what state vector components are going to be used in the model output. 
PD, (E) is the transition matrix of the noise process. 
U(t) is the white noise process vector. Includes the knowledge 
of f). 
h(E) is the noise process output. 
h(t) is the output vector of the noise process. 
The central problem of this Investigation can now be clearly defined. 
The information received by the radar and sequentially stored is the sum, 
Ан ДОГА (АР CER 
where y(t) represents the true distance and (7) tli« rior added to 
that. We will seek a procedure of machine computation to obtain yt) 
at each instant (-1,2,3,.... and, in a more general fashion, to obtain 
some χώ)- y(0,e, (O) where the error e, (t) is minimized in accordance 


with the mean-square error criterion. 


* ; Е 
Notice here that as we are treating the X-axis components this 
will have to be interpreted as distance along the X-axis. This apparent 
confusion could be avoided but we decided to maintain completely Magill's 


notation to provide some unity between his report and this application, 
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3. Consideration of the theory 


Later on we will study Kalman filter's adequacy to radar data filter- 
ing and will concentrate in the areas of application. But, obviously, 
a clear understanding of the theoretical resuits to be used has to be 
provided. This is the present chapter's objective. 

The elements required for a solution may be concentrated in the 
four equations presented in the last chapter which describe the message 
and noise processes. 

A simplification of this can be made and consists of joining to- 
gether the message and noise processes, reducing these to a single 
process which has as output the sum z(2)=- Y “Jn(() This is accom- 
plished by augmenting the state vector. 


According to this we now define: 


£) 
х@)= 1677 - state vector of the observable process. 


GE OLE о 


tea) transition matrix of the observable process. 


O Pwl 
2 (INA 
20 e 271,06 
vé)? ( PE) : 


γι 0/ 72,1117, 
md =(r I 2 - output vector of the observable process. 


Defining also, 


MUGAZ - covariance matrix of the inputs. 
ο ο τῶν» 
22) 2 Al) RE) 








The model for the process which generates the observed signal 
may now be reduced to the set of equations: 

Xx (t+) = é (Ü) x (É) + (tn) U(Ers) 
ZY = mx) , 
where U(t?) - U(t) » Of Fe/ 

Even if this group of definitions casts some shadow on the physical 
meaning of the parameters we will not try to clarify it right now. The 
next chapter will provide full understanding by means of an example. 

Ме may now represent the whole process in a block-diagram form, 
Fig. 1, where the parallel lines represent flow of vectors instead of 


a scalar value as in the line at right for z (7. 





vtr) Е | (9%) z0 


Fig. | Model of the observable process 


-2 
*As Qt) is a diagonal matrix, ? ΚΖ, is equivalent to apply 


the exponent -і to all non-zero elements of this matrix. 








Let us proceed now to look at the solution of the filtering problem. 

The objective is to find a recurrence formula to obtain for each 
sampling point the best estimate of the state vector, given both the 
previous one and observed state vectors at that sampling point. 

The best estimate of the state vector we need is by definition, 

X(t/t-1) 2 Ef x(t)/2,,}, 
Where Z,., represents all the available data up to and including time Z-/ 

(104) 2 z(9-(tr-), 
where z(é) is the observed signal at time Č and 2 (t/t-1) is the best 
estimate of the observed signal at time £ given the previous best 
estimate. It is possible to derive [3] a recurrence formula, 

X (E/t) =X (tht-1) είς) ἔ(ε/έ-) , ἐ12,5....., (1) 

where K(£ €) represents a variable gain vector. 

We now observe that the filter produces the best estimate at time C 
given the data collected up to and including that time by means of the 
following operations, 

a) Evaluates a best estimate at time Z given data collected up to 
and including Z-4 . This is simply done by multiplying the transition 
matrix 2(2-/) Бу the best estimate X (-//#-/) Р 

b) Corrects x (¢/t-1) by the addition to it of a term containing 
the information received at C . This term is the product of a weighting 
coefficient - the gain k(t, Е) - and the error in estimating the observed 
signal given data collected up to and including Z-/ . 

From what has been said we understand that the filter is essen- 
tially formed with a model of the message process and the added feature 
of being able to correct the updating of its output with the information 


received. 








The recurrence formula (1) is the solution to the filtering problem. 
Of course the means of updating the gain matrix is still not uncovered. 
The derivation of a recurrence formula for this matrix starts out 
considering that the best estimate can be obtained by minimizing the 
trace of the following matrix, 
РЕ)? Е {Хх (üt) ο 
where x (tlt) 2 X (2) -X (t/t) , 
X (t) being the true state vector of the message and noise gener- 
ated. 
The matrix P(t/t) is the covariance matrix of the estimation error 
and, obviously, gives an indication of the error value at any Z. 
Further defining, 
r? (tlt)? mY P(tlt-+) ro (9= ver [Z(E/E-Y (2) 


it is possible to obtain 


k (t,t) РИ m (2 
(2,2) TEE) (3) 


A recurrence formula for Р(2/2-) must still be obtained. Again 
Magill derives it based on the assumptions previously made and reaches 
the result, 

— le NT -ACE t)rr'le)]ele/t-)PT РИ) @ 

The recurrent procedure ¡s now completely obtainable, provided the 
initial covariance matrix РЕ previous knowledge) i s given. This is 
done by sequentially evaluating (ZZ), Х(Ё/0), 0, ЈА) , etc.. 

With the equations we have introduced, a filter block-diagram is 


possible. 








z(t) ; з — (Е x= X ft 





Fig. 2 Model of the optimal filter 
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h. Applications 


In this chapter we intend to analyze more closely the assumptions 
made before getting the results of the last chapter and test them in 
the problem of filtering radar data. We will then decide on the 
applicability to our problem of the theoretical results presented before. 

As it was sald, all Magill's derivations are based on the fact that 
the processes Involved are either Gauss-Markov or can be obtained from 
these by application of a linear matrix multiplier. Can the state 
vector of any flying target, whatever dynamic properties it may have, 
be considered a Gauss-Markov process? In fact, given any instantaneous 
state vector, it is always possible to consider some other previous 
time when the definition applies, i.e., 

[f CIC, x 92), J= pP ӘУ (е), lk, 
where Á is finite. 

In a physical sense we may say that if we want, at any instant Û, 
to obtain the state vector at É#/ , we may discard the knowledge we have 
of the target's state vector before the instant WAK. For example, if 
we consider an aircraft at £ , clearly it is seen that the knowledge 
of all its positions since the flight started is normally irrelevant to 
determine the next one but we are almost certainly interested in the 
last two or three. Of course the number K will vary with the type of 
target considered but in any case is a finite number. 

In Chapter two we said in d) that we would need a dynamic model 
for the tracked target. This model and the information of e) (some 
statistics of the inputs) and f) (statistics of the radar error) were 


all considered in the formulas then presented. 








Let's analyze first the message process, 
s (t+1) = 6 (t)s (t)+Q (tu (t+) 
y(t) - n" (9s (2) 
If we consider in the state vector only distance, speed and 
accelleration and this way eliminate the possibility of a third deri- 


vative in the motion we have, 


Ξ(ε | dl 6, | |501 |8,69 2,699,692 |6,69 
Ss |&,( 4,09 Bn l| |R| + 2,069 2 6) 5. 9| |U e) 
AED] US," AV) (SO) 92,4% 9,62 2,9 | 6,0) 

Simplifying the notation and multiplying matrices we obtain, 

5, = #, S, * ó, S; * #, 55, 2,0, + 2, 00а + Dsn бу 

The first three terms of the second member of the equation repre- 
sent the formula to obtain the distance at some instant Zv/ , the dis- 
tance, velocity and acceleration at £ being known. The components Φ, ; 
2, and 9, could be easily gotten through the Taylor series expansion 
and as we know they are 

@,=/, δι.” ὃν - 27”, 
where 7 is obviously the sampling period. 

The other components of the Ø matrix are easily obtained in the 
same fashion. 

Ме can now give full significance to A (tJ (tu) if we recall 
the definition of the U (tri) matrix given before. It is simply the 
way of correcting the next value of the target's state vector by adding 
the possible random inputs. For example, we may obtain the target's 
next velocity value if we know the present one, the acceleration and 
the sampling period but notice that we still have to consider the intro- 
duction, during the sampling period, of some wind velocity and then we 


have to complete the value of S Gel)= $, (в) + r s, (O by the addition 
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of U, (tr) 


The random inputs U, (t) and "NÉ, account for instantaneous 


changes of position and for variations of acceleration introduced by 
differences in the engine's thrust. 


All put together in the matrix equation we have, 
S, (t+) I T 4r? 5,(ἀ) | |1 
5, (ём) о / 


o o| Vs Ce) 
А s)| lo 
S») оо 7 


го | |0 (1) 


sS, 10 ох (е) 
One feature of the random inputs when considered together is that 


one of them may substitute the entire vector. Consider, for example, 
Uss (м) . Through this acceleration input we can account for any 


variation in position or velocity which instantaneously may happen. 


Looking now at the second equation of the message process and 


considering y(t) as a vector containing only distance, 


[r@] = o 


ο 


In the noise process we have a different case. We consider the 
noise at any RTT time to be independent from any other noise value. 

This makes all the elements in 2, (t) null. 

Let us make it clear also that the state vector W (t) has only 
one component which is, obviously, the error introduced in the measure- 


ment of distance, the only one made. The input Uy (tr) is then the 
actual error W(t+/) . 


We have then, 


[vct] - [o] [wo] * [i] [ο (ey), 
еј гне 


(5) 


(6) 
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A point we still did not look at is the gaussian quality of the 
random inputs. Can we admit that the acceleration imparted to an 
aircraft when considered in discrete amounts constitutes a random 
variable with gaussian distribution and zero mean? If we think of 
all the aircraft and all its flights we, naturally, would be ready 
to accept that. And actually, this is what happens once the filtering 
must be appropriate for any aircraft that can appear in the radar scope. 
This has to be kept in mind when values are going to be given to the 
inputs. We wil?! return to this problem in the next chapter. 

Having formed the message and noise processes we may now proceed 


to construct the model for the observable process, as explained in 


Chapter 2. 
So, 
23 
Ж(Ё/= 2 
e Sy (8 
м (2) 
[| T фт? ο 
ο / 7 O 
20)= О O / O 
O O O / 
/ о O O 
O / 0 Ó 
26% O O / O 
O о O / 
y 


(0 

νε τς (6 
Dy, Ct) 

mi [/ o o /] 








and with simplified notation, 


v d 6ξε ya = 5 2 
U 
2 & “5а s> “Se $2 
QC z) a= (5, 4; d Us 5 Us tw 

Us, Uu Us, Um ott, Un 

As the random variables 4,4, , Сз, Πᾶνα zero mean, we 
see that, 
2 
ЖЕМЕ Var 6, Я 

апа, 


£ J, 4, } = СЦ, дг)=о 


since %,, 6. , etc. are uncorrelated. 


г ооо Grud = o © δ 
jw on  oNKO O (vor, )% ο. 

М | ο ο 1 ο O O (ery)? ο 
O o o / © O ο (όν) 
σι ο ο о 
Ô τω O 0 
© o бф, о i 
O O о (4, 


where ( represents the standard deviation. 

With the matrix just obtained for 22 we are able to use the 
recurrence formulas presented in the preceeding chapter. 

But we still need the initial covariance matrix without which we 
cannot start the computation. 

|t can be observed very easily that at г, , i.e., when the first 
measurement of the observable process is made, the correlations between 


distance and speed errors, distance and acceleration errors, etc., are 


all null. We have then, 


Var 

Pl. [Previous knowledge) = ο WX. О 0 
9 
Ó 
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Of course the diagonal values must be based on the information 
about the target which we have even before we begin tracking it. To 
avoid a purely conjectural evaluation of the variance we may simply 
discard the three first received positions from the radar, as far as 
filtering is concerned, and make use of them to estimate the velocity 
and acceleration of the target. At Ù+ we will have then the position, 
the velocity and the acceleration of the target with errors the variance 
of which we can determine. 

Let us see what are those variances. Supposing that we know the 
variance of the error in the distance measurement: 

By definition Х,(&+2)= х,( 42)— % (tere) 

Then Var X, (Te+2) = Var [x (2.42) = x, (t.+2)] 

if we take X, (6 2) as the observed value, then we have, 

Ur X, (kr2) = r Y 
Ww 
because the error in the distance measurement is the error introduced 
by the radar. 
Ме want now Var Х, (622) . Using the last two points we have, 


using the same assumptions as above, 


xa (kra) - Xoltor2)- X tor) X (O2) (тан) 
— TA P 
e X r2) - 2) Xx (su)-X, GS) . 
" x 7 
a Х, (+2) — X, ( Zer/) 


7 
Finally Vor Xe (lore) < 2. var [X 6% (6807 = 
= = Vor 6, 


Let us obtain now 


Var X, (6+2) 








The three points we have permit us to write, 


Xs (+2) = X3(b,42)- Xz (ter2) = 
= _X2 (s + 22 = Xa (Ax) = Xs ( —— = x (t+) E 
7 7 
хк) - X, (hw). Xen) X). Xen - Xue) e (Er) X) 
A A A л μμ 


7 Z 
— x Gr2)- 2906 ы) _ X, (рә) -2X, (bv) +X; (25) 
7 2 y 2 
Now Ж х. (Lot2Z) = L Гин aA (br) ver X, CE 
= е, Yer Uz, 








5. Computation of the filtering process 


In the beginning of this work we showed our interest in the 
eventual machine computation of the recurrent process. It is now 
time to look at the procedure introduced in the last chapters and try 
to evaluate how effective the machine computation can be made. 

Im accordance with what we mentioned in chapter three the 
formulas below constitute the recurrent process, 

P(E)t-) » (t) ] T- eu t) e-ylee-) 8 C=)» ТН) 

ПОЛЫ) = m (1) РЕЛ) (Е) 


K (Ut) « FE) meg 
c *⁄¿ 7-7 


ВО) m™ (9 (t-Y (ta / ev) 
X (U) = BUDS (Ct) 2 уы (6 2) 


The program for this computation can be done very easily. We 
have in the next page a very simple flow-chart that illustrates 


this array of computations, 








START - 


yerror(2)= | 2(/-)- 






Fig. 3 Flow-chart for recurrence formulas 
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The problem of how to start the recurrent process has still to 
be considered. We said before that we could obtain the value of the 
initial covariance matrix by simply taking the three first points of 
the raw track, evaluating with them the position, velocity and accel- 
eration at the third point and, at the same time, calculate the initial 


covariance matrix of the state vector thus obtained. 


The formulas then mentioned are repeated here for convenience, 


Х, (6+2) =Уь»2) 
Х. (52). У =2)- yG@s*2 
7 


Y (Co +2) -Mi 1) _ Xll (У a) 
Xy A — πε ΓΝ. 


+ 
-- 


X,4(br2) = О 
Var X, (Lot2) « Var Usw 
Var X. (619 « 2, Yor Um 


И» Xz (612) - £ νέν» Um 


Vor X, (6:2J£ Var б» 


Considering the computations above we obtain the flow diagram 


on the next page. 


20 








X(3)- Y(3) 


С Formobe derived 








P// (3) = Q44 
P22[3) E or... 
P33 (Ze °° 
P44 (3) = 244 






Fer olas een) ved 





Fig. à Flow-chart for the optimum filter 
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This way we use the initial covariance matrix at Cot2 (Z= z) 
to obtain the next one at Z,*3 and use the observed positions at L +2 
to obtain a best estimate at Z*3 . Our filter will then begin filtering 
at and including 4*3 

To avoid too long statements which would add only to confusion in 
the program we divided the computation of the covariance matrix in 
parts as It can easily be seen. 

Another point in need of consideration is the way the program is 
implemented without the help of sub-routines, even if some matrix 
operations are repeated several times. We thought that some computing 
time would be saved by eliminating the call of sub-routines. This way 
there is no great compilation time spent and probably the compilation 


done comes out more efficient. 
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6. Example of a simplified simulation 


In this chapter we are finally ready to implement the techniques 
which have been studied until now. Ideally that should be done in a 
real case. A seemingly feasible way would be to obtain radar data from 
maneuvering aircrafts, project the points in the X,Y and Z axis and 
try to filter the three tracks thus obtained. This way of implementing 
our methods is, however, not too flexible. 

A better way is to create a track, add some noise and try to 
filter this out afterwards. This approach will provide us with the 
possibility of varying the track at our wil! and thus permit the 
evaluation of the filter's suitability in various conditions. Ном 
much the results can be relied upon depends, obviously, on the simu- 
lation's fidelity. This we will try to assess in the next chapter. 

As we have said before, to study our filtering problem it is 
irrelevant to consider more than the projection of the target's track 
on one of the coordinate axis. The procedure for the other axis is 
exactly the same. 

Limiting ourselves to the projection on the X-axis we can have 
three different situations - constant velocity, constant acceleration 
and sinusoidally varying distance. These correspond in three dimensions 
to constant velocity, constant acceleration and circular track. Any 
other track of irregular form can be reduced to a group of circles with 
accuracy as high as we want. 

It seems then only reasonable that we study tracks where these 
three possibilities happen in different orders and with different 
values. From here we may eventually conclude that the filter provides 


acceptable results fm any conceivable simulated situation. 
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Furthermore it seems logically acceptable to conclude the following: 

a) The filter will yield good results with any constant velocity 
given if that is the case for the lowest and highest expected velocities 
simulated. 

b) The filter will work in good conditions for any constant 
acceleration given that it works for the highest possible and fulfills 
the requisite a). 

c) The filter will definitely accept any simulation and yield good 
results if, furthermore, it gives acceptable results with the tightest 
circular maneuver possible = іп the sense that the acceleration on the 
X-axis reaches the highest value possible with the speeds available. 

At this point we will have to decide what types of targets are 
we going to be handling. Let us confine ourselves to airplanes. 

To simulate a track we will suppose an aircraft going-through all 
the situations presented above. 

Furthermore it seems desirable to have all the combinations of 
those pieces of track that seem to introduce different degrees of 
difficulty. 

The tightest maneuver we want to explore is the one which gives 
the greatest maximum acceleration of the motion's projection on the 
X-axis. It is easy to see by Fig. 2-29 in one of the references [4] 
that the highest speed and greatest bank angle will give the largest 
value for acceleration. Considering that the maximum possible man- 
euvering speed for the aircraft is 1000 knots and the banking angle 
70 degrees, we have that the projection on the X-axis of a circle 
described in the X-Y plane is, 


X= K Cos (Rot x t) 


where Ris the radius of turn in miles, Kot the rate of turn in degrees 
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per second. 


Extrapolating in the tables for 1000 mi we have, 


R = 32,000 Ug. οσο δις ο 


6,000 


Rolls 5.6 κ P = .05237 red /sec 


Then the maximum acceleration fs, 

e = 2: (pot)? = 6-6x(-05237)° ri] sec? = 
e € x- 0027€ mif sec? a «QUAE ж», /ғес? 

Let us assume that in a straight line the maximum acceleration 
possible is 1.805 g's. This gives, 

1,805 X32.2 J sec I Mile? = 983 x10 rm есе 

To have the largest range of speeds we wil! consider that our 
aircraft can fly from 10 to 1000 mi/hr, which is unreal but will give 
greater scope to our results. 

Let us see now what types of tracks will be worthwhile to exper- 
iment with. Below are the different combinations which seem to present 
different degrees of difficulty to our filter. 

a) Lowest speed, acceleration to highest speed, highest speed, 
circle at highest speed. 

b) Lowest speed, acceleration to highest speed, circle at highest 
speed, highest speed. 

c) Inverse sequence of a). 

d) Inverse sequence of b). 

Another factor that may be important is the time during which the 
aircraft will stay at the lowest and highest speeds. 

Two problems are still not solved. Опе îs the addition of noise 
to the track we are going to create. The other is the sampling interval 


which corresponds to the radar rotation period. 
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The noise will be considered to have a Gaussian distribution with 
a variance of four miles constant for any distance. Of course this is 
a rather strong simplification and in reality the variance will depend 
very much on the distance. We do not think that it is too important 
to deal with this problem right now. |n the next chapter this will 
be discussed again. 

The sampling period will vary between some limits in a real 
situation. But, again, that fact does not seem to be important in 
this simulation. 

The programming of the radar data simulation is very clear and 
can be analyzed in the Appendix. 

First, a time array is generated, МЕ (К). Then, for each piece 
of track we proceed as follows: 

For each sampling point the noise is generated with the RNDEV 
sub-routine which provides us with random numbers obeying to a 
gaussian distribution with variance of unity. DEV is the output 
number, which we multiply by МСОЕ to get a variance of valve COE. 
Then we evaluate the signal output vector using formulas of elementary 
physics. Finally, the distance is added to the noise giving the input 
to our filter. 

A last question is what values to give to the Q matrix. First 
of all we are going to assume that all the random inputs are null 
except the acceleration. Then we will try, by the use of several values 
for Q33 , to obtain one which fits best the whole filtering for all 
the cases we are going to experiment with. 


The values that finally were reached are printed out in the Appendix, 
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To help in the evaluation of the filtering quality a smoothness 


index was defined as, 


smoolh _ Var (ANT) —Var(K1E) 9 
— — x /00 4 


27 








7. Conclusions and acknowledgements 


The results obtained indicate that the Kalman filter cam be used 
in the filtering of radar data with some success. 

The cases taken as examples did indicate that in a wide range of 
real situations the filter will provide useful filtering, as far as 
the simplifications made are not the reason why the filter worked. 

The constant sampling period which was used could have been varied 
along the track without any difficulty. To prove this we only note that 
for each computation made after the arrival of a new radar position the 
value used for Z îs not related in any way to the last value of 7 . 

The radar error may have a bias depending on the distance. This 
bias has to be experimentally determined and its value subtracted from 
the range measured. If not it will eventually introduce a bias in the 
filtered values. 

The Gaussian quality of the radar error is a point of doubt. Only 
in the presence of a definite radar system can that be asserted. Ном- 
ever in a general way we observe that al] components of the radar noise 
seem to have a Gaussian quality. 

The value we should give to the Q matrix is critical. Ideally, 
for each type of target (helicopters, airplanes, missiles, ships, etc.) 
we should have a Q matrix. The values for Q33 can be chosen experi- 
mentally. The filter will work with less error when the variance of 
random inputs is smal], as it would be the case of a missile when come 
pared with an aircraft, especially if we consider a ballistic missile. 

The components of the motion on the three coordinate axis can 
obviously be evaluated in succession and with a normal speed computer 


the computation wil] be completed before the next set of data arrives. 
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X-SCALE - 1.000402 UNITS/INCH, 
Y-SCRLE = 1.007401 UNITS^INCH. 
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“SCALE = 1.00£2+92 UNITS INCH. 
Y-SCRLE = 1.00£+89 UNITS“ INCH. 
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X SCALE = 1.000202 LINITS^INCH. 
Y-SCALE = 2.007409 LINITS^INCH. 
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A SCALE = 1,002+02 UNITS INCH, 
Y SCALE - 2.00£+01 UNITS INCH, 
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XSCALE = 1.88E+82 UNITS- RYCH. 
Y-SCRLE = 1.00Е+01 UNITS- INCH, 
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А-ЭСАГЕ = 1.00E+02 UNITS INCH, 
Y-SCALE = 1,00£+08 UNITS/INCH, 
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K-SCALE = 1,082-+02 UNITS-TMNCH, 
‘(SCALE = 2.0024+41 UNITS INCH, 
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X-SCRLE = 1.00E4+02 UNITS INCH, 
Y-SCALE - 1.00E401 LINITS^INCH. 
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K-SCALE = 1.08£+82 UNITS/INCH, 
Y-SCALE = 1.007209 UNITS- INCH, 
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K-SCALE = 1,00£+02 UNTTS- INCH, 
Y-SCRLE - 2.09E+01 UNITS INCH, 
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K-SCALE = 1.88E+82 UNITS INCH, 
f-SCRLE = 1.007401 LINITS^INCH. 
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K-SCALE = 1,08£+82 UNITS/INCH. 
*-SCRALE - 1.007409 LNITS/^/INCH. 
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-SCALE = 2.00£+90 UNITS, INCH, 
CARLOS P. SIMOES MS THESIS 
OPTIMAL FILTERING OF RADAR DATA 
























Case c) 









600 991 092 083 994 


K-SCALE = 1.08£+82 UNITS INCH, 
Y “SCALE = 2.00£+91 UNITS INCH, 
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